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Abstract 

We prove the global existence and uniqueness of strong solutions for a 
compressible multifluid described by the barotropic Navier-Stokes equa- 
tions in dim = 1. The result holds when the diffusion coefficient depends 
on the pressure, ft relies on a global control in time of the L 2 norm of the 
space derivative of the density, via a new kind of entropy. 

§1 Introduction 

In this paper we show the well-posedness of a global strong solution to a mul- 
tifluid problem over R + x R characterized by the one-dimensional compressible 
barotropic Navier-Stokes equations. That is, we consider the following system 
of equations, 

d tP + d x ( P u) = o, (i.i) 

d t {pu) + d x (pu 2 ) + d x p(p,p) - d x {v(p,p)d x u) = 0, (1.2) 
d t {pn) + d x (pup) = 0, (1.3) 

with initial conditions given by: 

P\t=o = Po > 0, pu| t=0 = to , M|t=o = Mo- 

The conservation of mass (1.1), conservation of momentum (1.2) and conser- 
vation of species (1.3) describe the flow of a barotropic compressible viscous 
fluid defined for (t, x) £ R + x R. Here the density is given as p, the veloc- 
ity as u, the momentum as m, and the mass fraction p denotes the relative 
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weighting for each fluid component of the adiabatic exponent 7(/u) G R associ- 
ated to the generalized pressure p(p, fi) , thus effectively tracking the "mixing" 
of the fluid components. In one dimension the shear viscosity and the bulk 
viscosity collapse into a single coefficient function depending on p and p which 
we denote here v(p,p). Monofluid one-dimensional compressible Navier-Stokes 
equations have been studied by many authors when the viscosity coefficient v 
is a positive constant. The existence of weak solutions was first established by 

A. Kazhikhov and V. Shclukhin [55] for smooth enough data close to the equilib- 
rium (bounded away from zero). The case of discontinuous data (still bounded 
away from zero) was addressed by V. Shclukhin [86, 87] and then by D. Serre 
[85, 84] and D. Hoff [47]. First results concerning vanishing intial density were 
also obtained by V. Shclukhin [88]. In [49], D. Hoff proved the existence of 
global weak solutions with large discontinuous initial data, possibly having dif- 
ferent limits at x = ±oo. He proved moreover that the constructed solutions 
have strictly positive densities (vacuum states cannot form in finite time). In 
dimension greater than two, similar results were obtained by A. Matsumura 
and T. Nishida [65, 66, 67] for smooth data and D. Hoff [48] for discontinuous 
data close to the equilibrium. The first global existence result for initial density 
that are allowed to vanish was due to P.-L. Lions (see [61]). The result was 
later improved by B. Desjardins ([21]) and E. Feireisl et al. ([14, 34, 35, 36]). 
The class of solutions was then extended by A. Zlotnik, G.-Q. Chen, D. Hoff, 

B. Ducomet, and K. Trivisa in [15, 17, 28, 99, 107, 108] to the case of a thermally 
active compressible flow coupled by the systems chemical kinetics, where global 
existence results are shown for an Arrhenius type biphasic combustion reaction 
tracking only the reactants level of consumption. Y. Amirat and V. Shclukin 
have further provided [1] weak solutions for the case of a miscible flow in porous 
media. 

The problem of regularity and uniqueness of solutions was first analyzed by 
V. Solonnikov [89] for smooth initial data and for small time. However, the 
regularity may blow-up as the solution gets close to vacuum. This leads to 
another interesting question of whether vacuum may arise in finite time. D. 
Hoff and J. Smoller ([50]) show that any weak solution of the Navier-Stokes 
equations in one space dimension do not exhibit vacuum states, provided that 
no vacuum states are present initially. 

Interfacial multicomponcnt flows have been extensively studied in the litera- 
ture, and span a rich array of applied topics with natural analogues in continuum 
dynamics. For example, there has been numerous work on multicomponcnt flows 
in biological systems, including bifurcating vascular flows [19, 39] and pulsatile 
hemodynamics [54], in vitro tissue growth [59] and the "amoeboid motion" of 
cells by way of surface polymerization [31], chemotactic transport [37] in aque- 
ous media under chemical mixing (e.g. varying relative concentrations) applied 
to specialized cell types [30], as well as biological membrane dynamics due to 
local gradients in surface tension caused by flux in local boundary densities 
[83]. Another important and popular field of application is that of dispersed 
nanoparticlcs in colloidal media (e.g. aerosols, emulsifications, sols, foams, etc.) 
[68, 70, 98] applied to, for example, electrospray techniques in designing solar 
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cells [41], or more generally to diagnostic and flow analysis in the applied mate- 
rial sciences [11, 57]. In addition, phase separation and spinodal decomposition 
have received a great deal of attention [4, 13, 62], especially with respect to 
morphological engineering [45]. Another field which is heavily weighted with 
multifluid applications is that of combustion dynamics [102] and chemical ki- 
netics [42, 52], where the conservation of species equation (1.3) is regularly 
invoked, including numerous topics in reaction diffusion dynamics and phase 
mixing, spanning many essential topics in the atmospheric [20, 44] and geo- 
physical [81, 101] sciences. In electrochemistry and chemical engineering recent 
work has been done on porous multiphase fuel cells [82, 80], and in sonochcm- 
istry recent studies have shown acoustically induced transport properties across 
intcrfacial phase changes [25, 73]. Finally in the fields of astronomy and as- 
trophysics exotic multicomponent magnetohydrodynamic plasmas are studied 
[72, 78]. 

Many applied results exist for computational methods and schemes for solv- 
ing multicomponent flows. Let us briefly mention some notable examples. An 
early generalized numerical approach in multiphase modeling was presented by 
F. Harlow and A. Amsden in [43], which provides an extensive system of dynam- 
ically coupled phases using a conservation of species (1.3) equation obeying a 
number of relevant physical boundary conditions and which applies, in particu- 
lar, to compressible flows. In [29] J. Dukowicz implements a particle-fluid model 
for incompressible sprays, an approach extended by G. Faeth in [32, 33] to com- 
bustion flows. D. Youngs then, in [106], expanded numerical mixing regimes to 
include interfacial turbulent effects. These basic schemes and approaches have 
been applied by a large number of authors to a large number of fields, modeling 
an extremely diverse number of natural phenomenon, from star formation [53] 
to volcanic eruptions [77] . Some good reviews of the foundational numerics of 
these approaches can be found in the books of C. Hirsch [46], P. Shih-I and L. 
Shijun [79] and M. Feistauer, J. Felcman, and I. Straskraba [38]. 

Let us briefly outline the physical meaning of the subject of this paper, 
namely, the system of equations (1.1)-(1.3). Here we have a barotropic system 
with the flow driven by a pressure p that depends on the density p and the mass 
fraction p of each chemical/phase component of the system. Since the function 
7(/x) depends on the constant heat capacity ratios 7$ > 1 of each component of 
the multifluid, the pressure p(p, p) effectively traces the thermodynamic "sig- 
nature" of mixing chemicals/phases in solution. Note that this is very similar, 
for example, to the system of equations set out in [102], except here, for sim- 
plicity, we have neglected the associated diffusion and chemical kinetics which 
break the strict (and mathematically convenient) conservation in the species 
equation (1.3). Another important facet of the system (1.1)-(1.3) is that the 
viscosity v is a function of the pressure p. Much recent work has been done 
by M. Franta, M. Bulicek, J. Malek, and K. Rajagopal in providing results on 
these type of viscosity laws [12, 40, 63, 64]. Moreover, since the form of the pres- 
sure p is chosen up to any state equation that satisfies the assumptions given 
in §2, the formulation is general enough to include, for example, multi-nuclear 
regimes. That is, in addition to describing the flow of mixing fluids charac- 
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terized by their concentrations with respect to their heat capacity ratios, this 
construction also educes applications in nuclear hydrodynamics, where one can 
derive the pressure law using the time-dependent Hartee-Fock approximation 
[103]. Such a nuclear fluid obeying the assumptions given by the Eddington 
Standard Model for stellar phenomena has a pressure law [34] that takes the 
form, p(p) = Cip 3 — C2P 2 + Czp 7 ^ where Ci, C2 and C3 are positive constants. 
In particular, this exotic pressure law can be shown to model nontrivial physical 
phenomena; such as spin and isospin wavefront propogation in nuclear fluids. 
It has further been shown to be in good agreement with nuclear hydrodynamic 
models of the sun [34]. Thus the result in §2 allows us to extend the above to 
nuclear multifluids that satisfies 

P {p,p) = Cip»V> - C 2 p»M + C 3 p^\ 

as long as it verifies the conditions given in §2. It however remains to be seen if 
quantum multi-molecular fluids [104] have an analogous formulation. 

At the level of the mathematical results, incompressible (for general back- 
ground on the incompressible Navier-Stokes equations see [60]) multicomponent 
flows have been addressed by a number of authors. First S. Antontsev and 
A. Kazhikhov in [2], A. Kazhikhov in [56], S. Antontsev, A. Kazhikhov and 
V. Monakhov in [3], and B. Desjardins in [21] show results for mixing flows 
where homogenization of the density p is allowed. These solutions can be seen 
in contrast with P.L. Lions' and R. DiPerna's solutions in [22, 74] which pro- 
vide a multiphase solution for immiscible inhomogcnizable flows given discrete 
constant densities for each component. A. Nouri, F. Poupaud and Y. Demay 
[74, 75] extend these results to functional densities where boundary components 
dfli are set between each fluid domain fli that satisfy the so-called kinematic 
condition, which restricts the viscosity v to obey d t i> + u-\Iv = (see [76] for fur- 
ther discussion on the kinematic condition). These results apply to immiscible 
flows with boundary surfaces that effectively fix the number of fluid particles on 
the interface. These solutions were then further extended by N. Tanaka [96, 97], 
V. Solonnikov and A. Tani [90, 94, 91, 92, 93] to include boundary conditions 
tracking both the surface tension at the interface using a mean curvature flow 
on the interfacial surface, as well as the inclusion of self-gravitating parcels. 

In this paper we consider viscosity coefficients depending on the pressure sat- 
isfying a barotropic-type pressure law, a result based upon the paper of A. Mel- 
let and A. Vasscur in [69] and extended to the multifluid case with a viscosity 
functional v{p) given no a priori uniform bound from below. Thus, in addi- 
tion to modeling the miscible multiflow regimes that have generated substantial 
physical interest (see above), our result further incorporates a very inclusive 
form of the generalized viscosity. We show the global existence with uniqueness 
result for a one-dimensional compressible barotropic multicomponent Navier- 
Stokes problem. In order to aquire the existence result, we rely heavily on an 
energy inequality provided by D. Brcsch and B. Desjardins (see for example 
[6] and [10]). This beautiful and powerful tool is central to our result, and, 
as it turns out, the breakdown of this calculation is the only (known) obstruc- 
tion to acquiring similar results in dimension greater than one. Next we obtain 
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the uniqueness result by adapting a proof of Solonnikov's [89] to the case of the 
barotropic system (1.1)-(1.2) coupled to the species conservation equation (1.3). 

Let us take this opportunity to discuss difficulties and related systems of 
equations in higher dimension. Again, the present result relies heavily on the 
calculation of an energy inequality (see §3) as provided by D. Bresch and B. Des- 
jardins [5, 7, 9]. However, in dimension n > 2, the derivation of this entropy 
inequality leads to an unnatural form of the viscosity coefficiant v(p,ii); which 
is to say, the calculation no longer demonstrates the type of symmetry which 
leads to the essential cancellation of singularities (for example see [49]) required 
in the calculation (see [71] for the monofluid case). 

We briefly recall some exciting results known for compressible fluids in higher 
dimension, and note that extending these to multifluid regimes introduces both 
beautiful and difficult mathematics, while also addressing very important and 
physically relevant questions in the applied fields. For example, a result of 
A. Valli and W. Zajaczkowski [100] shows global weak solutions to the mutlidi- 
mensional problem for a heat conducting fluid with inflow and outflow conditions 
on the boundary. In [95] A. Solonnikov and A. Tani offer a uniqueness proof 
for an isentropic compressible problem given a free boundary in the presence of 
surface tension. D. Hoff and E. Tsyganov next provide a very nice extension of 
the system to find weak solutions to the compressible magnetohydrodynamics 
regime in [51]. G. Chen and M. Kratka in [16] further show a free boundary 
result for a heat-conducting flow given spherically symmetric initial data and 
a constant viscosity coefficient in higher dimension; a result which is extended 
by E. Feireisl's work [34] under the notion of the variational solution for heat- 
conducting flows in multiple dimensions; though this result restricts the form 
of the equation of state. Further existence results are provided by B. Ducomet 
and E. Feircisl [26, 27] for gaseous stars and the compressible heat-conducting 
magnetohydrodynamic regime. Further, D. Donatelli and K. Trivisa [23, 24] 
have extended the existence results for the coupled chemical kinetics system 
mentioned above to higher dimensions. In [70] A. Mellet and A. Vasseur pro- 
vide global weak solutions for a compressible barotropic regime coupled to the 
Vlasov-Fokker-Planck equation, which characterizes the evolution of dispersed 
particles in compressible fluids, such as with spray phenomenon. Finally, im- 
portant results of D. Bresch and B. Desjardins, in a very recent paper [8], has 
worked to extend the existence results to a more general framework (using their 
energy inequality) for a viscous compressible heat-conducting fluid. 

We conclude by noting a number of important and interesting results related 
to vacuum solutions. That is, though in this work we are concerned with densi- 
ties that obey uniform bounds in R, a number of nice results exist for the case 
where over some open U C R, 



which is to say, solutions that incorporate vacuum states. For example, T. Yang 
and C. Zhu show in [105] global existence for a ID isentropic fluid connected 
continuiously to a vacuum state boundary with a density dependent viscosity. 
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Additionally, in dimension one, a recent result by C. Cho and H. Kim [18] 
provides unique strong local solutions to a viscous polytropic fluid, where they 
utilize a compatibility condition on the initial data. 



§2 Statement of Result 

Let us first state the hypothesis we make on the pressure and viscosity functional 
p(p,p) and v(p,p). First we assume that the pressure p(p,p) is an increasing 
function of the density p such that a.e., 

d pP (p,p)>0. (2.1) 

The viscosity coefficient v(p, p) is chosen such that it satisfies the following 
relation, 

"(P. M) = pd P P(p, P)tp'{p{P, p)), (2-2) 

where ip(p) is a function of the pressure restricted only by the form of its deriva- 
tive in p. 

We consider a multifluid for which the pressure functional does not change 
too much with respect to the fractional mass. Namely, Consider two 7 > 1 and 
7 > 1, where 7 < 7 < 7 up to the constraint that, 

^ 1/2 <^, (2.3) 



7 7 

l ' 2 > i±m _ 1. ( 2.4) 



7 7 

These relations are satisfied when 7 = 1.4 and 7 = 1.3, for example. 

Then, we ascribe the existence of constants C > such that the following 
conditions hold: 

ip'(p) > Csup(p- a ,p- 5 ), 
P^/C < p{p, p) < Cp^ for p > 1, p e K, (2.5) 
P 1 IC < p{p, p) < Cp^ for p < 1, p e M, 

where a and a are such that 

X-^< a <X±i^, (2.6) 

1Z^> S >1±^-1. (2.7) 

7 7 

Note that the existence of a and a comes from (2.3) and (2.4). 

Next we set conditions on the derivatives of the pressure in p and p; given 
first in p by, 

p^- 1 /C<d p p(p,p)<Cp^- 1 for p>l,pGR, 
p<- x IC <d p p{p,p)<Cp x <- 1 for p<l,peR, 
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and in p by, 

9^(p,M)<<V for p>l, M eR, 
d^p(p, p) < Cp 1 for p < 1, p £ R. 

Notice that a simple pressure which satisfies these conditions is, p(p, p) = 
C(p)p J ^ where 1/C < C(p) < C and with two constants 71 and 72 such 
that, 

7 < 7i < 7(m) < 72 < 7- 

In particular note that (2.1), (2.2), (2.5), and (2.8)-(2.9) are quite general as- 
sumptions, while the strong conditions, (2.3) and (2.4), have the effect of con- 
straining the amount p(p, p) can change with respect to p. 

For the sake of clarity we define H 1 (R) as the space consisting of all functions 
p for which, 

/ (d xP ) 2 dx < C. 

JR 

This paper is dedicated to the proof of the following theorem. 

Theorem. Assume a pressure p(p,p) and viscosity v{p,p) satisfying (2.2), 
(2.5), and (2.8)-(2.9) where the adiabatic limits 7 and 7 verify (2.3)-(2.4), and 
take initial data (p ,u ,p ) for which there exists positive constants ~g(0) and 
g(0) such that 

< g(0) <p < £»(0) < 00, 
Po € H 1 (R), uq £H 1 (R), p £H 1 (R), 

S"(p ,Po)dx < +00, 
\d x p \ < Cp , 

where § is the internal energy as defined in (3.2). We additionally assume the 
existence of constants R,S>0 and p, p > where po = p for \x\ > R and 
po = p for \x\ > S. Then there exists a global strong solution to (l.l)-(l.S) on 
M+ x M such that for every T > we have 

perfO.Tifl'fR)), d tP £ L 2 ((0,T) x R), 
u£ i oc (0,T;iJ 1 (R))n J L 2 (0,T;iJ 2 (R)), d t u£L 2 ((0,T) x R), 
p x £ L°°(0, T; L°°(R)), d t p £ L°°{0, T; L 2 (R)). 

Furthermore, there exist positive constants g(T) and ~g(T) depending only on T, 
such that 

< g(T) < p(t, x) < g(T) < 00, V(i, x) £ (0, T) x R. 

Additionally, when ip"(p), d pp p(p,p), and d pfl p(p, p) are each locally bounded 
then this solution is unique in the class of weak solutions satisfying the entropy 
inequalities of §3. 
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It is worth remarking that our results are slightly stronger than those pre- 
sented in the statement of the theorem above. Namely, the conditions po = p 
for \x\ > R and p = A f° r \ x \ > S with respect to constants R and S can 
be relaxed, such that simply choosing po and po close to the reference values p 
and p, is permissible as long as the internal energy £ (p, p) remains integrable 
at t = 0. 

We additionally use the existence of the short-time solution to the system 
(1.1)-(1.3), which follows from [89]. That is, as we show explicitly in §4, applying 
(2.5) and (2.8) to (2.2) provides that for every (p, p) the viscosity coefficient 
v(p, l 1 ) ^ C> for a positive constant C. This leads to the following proposition. 

Proposition 2.1 (Solonnikov). For initial data (po,Uo,Po) taken with respect 
to the positive constants g(0) and g(0) satisfying 

< q(0) <Po < g(0) < oo, 
PoeH^M), uoeH^R), poeH^R), 

and assuming that v{p,p) > C for a positive constant C, then there exists a 
T s > for each such that (1.1)- (1.3) has a unique solution (p,u,p) on (0,T S ) 
for each T r < T s satisfying 

p e L°°(0, T r ; H 1 (M)) 1 d tP e L 2 ((0, T r ) x K), 
u e L 2 (0,T r - 7 H 2 (R)), d t u e L 2 ((0,T r ) x R), 
G L^T^OR)), 9 tM e L°°((0,r r );L 2 (K)); 

and i/iere exists two positive constants, g > and g r < oo, sucft i/iai g < 
p(x, t) < g r for all t G (0, T s ). 

The proof of Solonnikov's proposition 2.1 as presented in [89] follows with 
the addition of equation (1.3) by applying Duhamel's principle to the transport 
equation in p given the regularity which we demonstrate in §4, in much the 
same way Duhamel's principle is applied to p in [89] for the continuity equation. 
The rest of the proof then pushes through directly by virtue of the calculation 
shown in §5 of this work. 

§3 Energy Inequalities 

In this section we derive two inequalities in order to gain enough control over 
(1.1)-(1.3) to prove the theorem. That is, from these inequalities we obtain a 
priori estimates that hold for smooth solutions and then prove the existence 
result using Solonnikov's short-time solution. The first inequality we show is 
the classical entropy inequality adapted to the context of a multifluid, while the 
second is an additional energy inequality derived using a technique discovered 
by D. Bresch and B. Desjardins that effectively fixes the form of the viscosity 
coefficient v(p,p). 



§3 Energy Inequalities 



9 



A simple calculation is required in order to obtain the classical entropy 
inequality (in the sense of [10] and [71]). That is, multiplying the momentum 
equation (1.2) by u and integrating we find, 

Here <o(p, p) is the internal energy functional effectively tempered by a fixed 
constant reference density p < oo and a fixed constant reference mass fraction 
A < C, given by 



, p) = P ( 
Jp 



p (p( s ,p) -p(p,p) 



ds+p(p,p) -p(p,p). (3.2) 



Let us make this calculation precise. First we restrict to the first two terms 
of (1.2) to notice that, 

d t {pu) + d x (pu 2 ) = pd t u + pud x u + u (d t p + d x (pu)) , 

where subsequently multiplying through by a factor of u and integrating gives, 

J ju 2 (dtp + d x (puj) + i (pd t u 2 + pud x u 2 ^j jefo. 

This can be easily rewritten using (1.1), as 

/ (udt(pu) + ud x (pu 2 )) dx = - j id t (pu 2 ) + d x (pu 3 ) Xdx. 
Jr * Jr J 

Likewise the pressure term p x from (1.2) is multiplied through by a factor of 
u and integrated. In particular, the form this term takes in (3.1) is derived from 
a pressure p(p, p) that satisfies a conservation law (shown in §4) for a tempered 
internal energy <§(p, p). To see this, first notice that for any function of p and 
p we have, 

/ d t £(p,p)dx = / d p S'(p,p)d t pdx+ / d li ^(p,p)d t pdx 
Jr Jr Jr 

= ud f _ t S'(p,p)d x pdx - / d p (S'(p,p)d x (pu)dx 
Jr Jr 

udpS{p,p)d x pdx - j pd p §(p,p)d x udx 
i Jr 

/ ud p S(p,p)d x pdx. 
Jr 
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But here, since d x S{p, p) = d p Sd x p + d^&dxp, we can write, 



R 



d t £(p,p)dx 




which gives 



Jr Jr 




(3.3) 



Using (3.2) we find 

£{p,V)= pd p S'(p,fi)+p(p,il)-p(p,p), where S , (p,p) = 0, (3.4) 
such that computing pd p S{p,p) — S^p.p) arrives with the desired equality, 



This internal energy <§ over K arises in [49] for the single component case, 
where there G(p, p') is set as the potential energy density and treated in a similar 
fashion. Note that as ${p, p) is tempered with respect to a reference density p 
and a reference fractional mass p, this is all that is needed to control the sign on 
the internal energy S{p,p), with the only qualification coming from §3 which 
gives cases on the limits of integration. 

It is further worth mentioning that the above terms comprise an entropy 
y{p, u, p) of the system (as well as the entropy term of inequality (3.1)), where 
we write the integrable function, 



The final step in recovering (3.1) is to calculate the remaining diffusion term, 
which follows directly upon integration by parts. That is, after multiplying 
through by u and integrating by parts we see that 




(3.5) 



y{p, u,p) = — + S"(p, p). 




which leads to the result; namely (3.1). 
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3.1 Additional Energy Inequality 

The following lemma provides the second energy inequality that we use in order 
to prove the theorem. 

Lemma 3.1. For solutions of (l.l)-(l.S) we have 

T+ I {^+ P - 1 d x ^(p)\ 2 + ^(p^)}dx+ f p-^'(p)(d x p(p, fJ ,)) 2 dx = 0, (3.6) 

providing the following constraint on the viscosity v{p,p): 

v{p,n) =pd p pip'(p). (3.7) 

Proof. Take the continuity equation and the transport equation in /i and mul- 
tiply through by derivatives of a function of the pressure ip{p) such that, 

d p ip(p){d t p + d x (pu)} = 0, 

d^(p)[d t p + ud x pj = 0, 

where adding the components together gives, 

d t tp(p) + ud x ip(p) + pd p ip(p)d x u = 0. 

A derivation in x provides that 

d t (d x ip(p)) + d x (ud x ip(p)) + d x (pd p ip(p)d x u) = 0, 

which we expand to 

d t (pp~ 1 d x ip(p)) + d x (pp^ 1 ud x tp(p)) +d x (pd p ip{p)d x u) =0, 

such that adding it back to the momentum equation (1.2) and applying condition 
(3.7) arrives with 

dt{p{u + p~ 1 d x ip(p)}) + d x (pu{u + p^ 1 d x ip(p)}) +d x p(p,p) = 0. 

Multiplying this by (w + p~ 1 d x ip(p)) then gives, 

1 2 1 2 

-d t {p\u + p~ 1 d x ip(p)\ } + -d x {pu\u + p^ 1 d x i}){p)\ } 

+ {u + p^ 1 d x il>(p)}d x p{p, p) = 0, 

which when integrated becomes 

1 / {^ + p- 1 d x ^( P )\ 2 + ^(p,p)}dx+ [ p- 1 <p'(p)(d x p(p,p)) 2 dx = 0, 
completing the proof. □ 



12 



Compressible Multifmid 



§4 Establishing the Existence Theorem 

In this section our aim is to apply the inequalities in §3 predicated on the 
formulation in §2 to acquire the existence half of the theorem. However, in 
order to do this we must first confirm that the energy inequalities satisfy the 
appropriate bounds. Let us demonstrate this principle for both (3.1) and (3.6) 
in the form of the following lemma. 

Lemma 4.1. For any solution (p,u,p.) of (l.l)-(l.S) verifying, 

I {po^T +£(po,(J<o)}dx < +oo (4.1) 
Js. 1 z ' 



and 



2 



Pol , d x 7p(p ) 



u + 



Po 



+ £(po,lJ,o)}dx < +oo, (4.2) 



we have that 



esssup / \p— + &{p, pL)\dx + / / v(p, p,)\d x u\ 2 dxdt < C, (4.3) 

[0,T] Js 2 J Jo Jl 



and 



esssup / (JL+^M 2 + g( p ^)\dx+ C ( t^El\d x p\ 2 dxdt<C. 
[o.T] Jm 2 1 p > Jo Jr P 

(4.4) 



Proof. It suffices if every term on the left side of both inequality (4.3) and (4.4) 
can be shown to be nonnegative. 

First notice that we clearly have that pu 2 > for any barotropic fluid over 
M, since p is strictly nonnegative. To check that §{p,p) > Owe simply refer 
to the definition given in (3.4). Indeed | p'^y^> j > when p > p and 

( gkfffifcgl ) < for p < p, which implies 



f. 

J p 



■p{*,n)-P(P,n) da>0m 



Together with (3.2) this gives that S{p,p) > . 

Next we check the viscosity coefficient v{p,p). Here the positivity follows 
from (2.2), where again the pressure is increasing in p satisfying (2.1) and the 
density is positive definite away from the vacuum solution (which we show is 
forbidden due to proposition 4.1), so for a ip'(p) satisfying (2.5) we see that 
i>'{p) > 0. Similarly, the last term on the right in (4.4) follows away from 
vacuum, where again we only rely upon the fact from §2 that ip'(p) > 0. □ 
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These results provide the estimates that we use for the remainder of the 
paper. That is, it is well-known (for example see Theorem 7.2 in [61] and the 
results in [69]) that the existence of a global strong solution to the system (1.1)- 

(1.2) follows by regularity analysis in tandem with (3.1) and (3.6). Below we 
present a similar approach for the case of a mixing multicomponent fluid (1.1)- 

(1.3) using only what we have found above; namely, that (3.1) and (3.6) provide 
the following a priori bounds: 

\\^IS(P, /J,)d x u\\ L 2( . T . L 2( R) ) < C, 

\\VPu\\l°°(0,T;L2(R)) < C, (4.5) 
IK(PjM)IU°°(0,T;L1(R)) < C, 



along with, 



\\(dx1p(p)/s/p)\\Loo( ,T;L2{R)) < C, 

\M(p)/p) 1/2 d xP (p^)\\LH n ,T^ m) <c. 



We will use these inequalities extensively for the remainder of the paper. 

As a remark, if we denote the internal energy density e as being characterized 
by the relations, 

f de(s,u) , . , . 

P e = P — a ds Wlth e{p,n) = 0, 

Jp dp (4.7) 

and d p e(p, p) = p~ 2 p(p, p) - p~ 2 p(p, p), 
then e is closely related to the specific internal energy e s , defined by 

_ f p p(s) 



-ds, 

s z 



which is provided for the single barotropic compressible fluid case in [76] and 
[34]; but in the multifluid context, since the internal energy S is tempered up to 
some constant reference density p, the usual form of the specific internal energy 
inherits a tempering in p as well, which is what is provided here by the function 
e. We also note that the tempered internal energy $ now satisfies the following 
conservation form as mentioned in §3, 

d t g{p,p) + d x {£{p,p)u) + |p 2 9 p e(p, p) + p(p, p) - p(p, p)}d x u = 0, (4.8) 

where it is easy to confirm that upon integration this recovers (3.5). 



4.1 Bounds on the Density 

For the existence theorem we need to establish a bound for the density in the 
space L°°(0, T; H 1 (R)). To achieve this we first establish uniform bounds on 
the density. 
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Proposition 4.1. For every T > there exist two distinct positive constants g 
and ~g such that 

8<p(t,x)<g V(i, x) e [0, T] x R. (4.9) 

Showing this proposition requires the following three lemmas which provide 
the groundwork for its subsequent proof. 

Lemma 4.2. Let & > be a function defined on [0, +oo) x R where JP(-,x) 
is uniformly continuous with respect to x e R and w/iere i/iere exisis a 5 > 
wi£/i J?(0,a;) > 5 /or any x. Then there exists an e > sucft £/iai for any 
constant C > i/iere exists a constant K > so that for any nonnegative 
function f verifying J R .*?(f(x),x)(lx < C, and for any x e R, there exists a 
point x\ G I = [xq — K, x + K] such that f{x\) > e. 

Proof. For any fixed & there exists a C such that for all y < e we have, 

ny,*)>± 

since J?(0,x) > S for any x and & is uniformly continuous in x. Let us fix 
C > and define 

K = 2CC. (4.10) 

We show that this K verifies the desired properties. Here we utilize a proof 
by contradiction in the spirit of [71]. Assume that we can find a nonnegative 
function / verifying J R J?(f(x), x)dx < C and an xq e R with 

ess sup / < e, 

KG/ 

where I = [xq — K, x + K]. Since & > this implies 

C"> / &{f(x),x)dx > [ &{f(x),x)dx> [ -Ldx, 
Jr J i J i 2C 

which yields C > K/C in contradiction to (4.10). □ 

Additionally we require the following technical lemma. 

Lemma 4.3. Providing (2.5) then (3.2) yields, 

p*< + Cp<C + l for p<l, 
pU^<C + C<?(p,Li) for p>l. 

Proof. Trivially, when p < 1 we have that p~' + Cp < C + 1. When p > 1 we 
use (2.5) to expand <S(p,p) where (3.2) gives that as p — > oo the p*> dominates 
such that scaling the constant correctly provides the result. □ 

Now we are able to find uniform positive bounds g and ~g on the density 
which inherently preclude the vacuum and concentration states. 
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Lemma 4.4. Assume that (2.3), (2.4) and (2.5)-(2.9) are satisfied and let 

dxZ{p) = M P <i}d x p- v + l {p >i } d xP a . (4.11) 

Then there exists an n > and a > such that for any K > there exists a 
Ck with 

R£(p)||loo(o,t ; li(i) < C K (4.12) 
for every loGt and I = [x — K, Xq + K] . 

Proof. First recall that the pressure satisfies 

d x p(p, p) = d p p(p, p)d x p + d^p(p, p)d x p. 

Here we are concerned with two cases, namely when p < 1 and when p > 1. For 
the case when p < 1 we multiply through by p~ 1 / 2 p(p, p)~— where a is given 
by (2.6), which yields 

dppjp, p)d x p = d x p{p, p) _ d^pjp, p)d x p 
y/PPfal*)- VPP(P^)- VPP(P^)- ' 

Likewise for p > 1 we multiply through by p(p, p)~ a p~ 1/>2 given a from (2.7) 
such that 

d p p(p, p)d x p = d x p{p, p) _ d^pjp, p)d x p 
^pp(p,p) w ^PP(p,p)" ^PP{p,^T ' 

In order to demonstrate the lemma we will control the right hand side of both 
(4.13) and (4.14) such that each is bounded in L°°(0, T; Lj oc (R)) for any point 
in M up to some fixed subinterval /. 

Towards this, we first show that p~ 1 d x p is bounded in L°°(0, T; L°°(R)). 
That is, take a derivation in x of (1.3) in order to write 

d t {pp~ 1 d x p) + dxipup^dxp) = 0, (4.15) 

such that multiplying through by a function ■d'(p~ 1 d x p) — i?' achieves 

dtip^p^d^)) + d^pw&ip-^dxp)) = o. 

Now choose such that for every test function d(p~ 1 d x p) € T>(R) with compact 
support, the function $ vanishes almost everywhere over the finite interval T = 
[—M,M], with M a constant. Upon integration this implies 



T d_ 

dt 



I pd{p 1 d x p)dxdt = 0, 
Jr 



such that for an appropriate choice of initial condition, where p Q 1 d x po G X, we 
find 



ess sup / pd(p 1 d x p)dx = 0. 
[o,T] Jr 
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This implies that pd(p~ 1 d x p) = almost everywhere for all (t, x) £ (0, T) x R, 
and so we can conclude that the argument of i? takes values over the interval 1, 
or more clearly that for p a.e. \p~ 1 d x p\ < M. This is then enough to educe the 
norm: 

||p _1 Qi:M||i,°°(o,t ; l°°(r)) < M - 

However, this is not yet enough to control the last term on the right for the 
two cases. In (4.13) applying (2.5) and (2.9) further provides 



^/pd^pjp, p) 
P(p, P)^ 



<C p 



7 — a 7 +7 



for 



for a positive constant Co- Using (2.6) from above we have that 7 > cry — 1/2, 
and so the positivity of the exponent gives 



C p7-«7+3 < C 



for 



P< 1, 



which leads to, 



{p<i}- 



y/pd^p, p) 



p(p,p)- 



< C. 



(4.16) 



L°°(0,T;L° 



Similarly for (4.14) we apply (2.5) and (2.9) to see that, 

J p • " ' 1 2 for p > 1. 



sfp~d»P{p,P) < ^ 7-«7+^ 



p(p,m) c 

Notice that since (2.7) provides 7 > 7 — S7 + 1/2, then applying lemma 4.3 
implies 

/ Vpd^pip, p) s 



4p>i> 



< C + C£(p,p). 



Integrating over / gives 



y/pdpP(p, P) 



{p - 1} PiP^r 

such that applying (4.5) establishes 

Vp~d»p{p, p) 



dx < 2KC + C [ £(p,p)dx, 
Jr 



{p>i}- 



P(P, Pf 



<c. 



K ■ 



L°°(0,T;L i 1 oc (R)) 



for Ck a constant depending only on K. 

Now consider the d x p term in (4.13) where here again we treat the two cases 
p < 1 and p > 1 separately. For the case p < 1 notice that we have by the 
bound on V>'(p) m (2-5) that 



{p<i} 



^P(P, p) 

Vpp{p, p)- 



c\i {p < 1} p- 1 / 2 d xP (p,p) 1 --\ < c|i {p <i } p- 1/2 a^(p)l- 
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Upon integration (4.6) gives 

/ \l {p < 1} p- 1 / 2 d xP (p,») 1 --\ 2 dx < C [ \l {p < 1} p-^ 2 d x tp(p)\ 2 dx 
Jr Jr 

<C, 

and so we obtain 

U{ P <i}P~ 1/2 p(p, p)~-d x p(p, m)I|l-(0,T;L=(R)) < C. 
Similarly for p > 1 we apply (2.5), giving 
d x p{p, P) 



li 



l ^^pp(p,»F 

such that integrating and utilizing (4.6) yields 



C|l { p>i } p- 1/2 0*p(p,AO 1_Q | < C\l {p > 1}P -^ 2 d x ^{p)\, 



< a 



and so 



l|l{p>i}P 1/2 P(P,P) a d x p(p,p)\\L~(o,T;mR)) <C. 



Combining these results we have thus acquired the important bound on the 
left sides of (4.13) and (4.14): 



l|l{p<i}P 1/2 P(P,P) -9 pP (p, P)dxP\\ L <*>(0,T;Ll oc 



+ I|1{p>i}P 1/2 P(P,P) a d P p{p,p)d x p\\ L ^ (S) , T . L i {I)) < C K - 



(4.17) 



It remains to show that for all p we have bounds on some power of the spatial 
derivative p x . First notice that when p < 1 applying (2.5) and (2.8) to the left 
of (4.13) provides 



d p p(p, p)d x p 



> qi { „< 1}P ™- 3 / 2 cM = C|l {p < 1}P 



{ ^ ir Vpp(p,P^ 
such that upon squaring and integrating wc find 



.7—07—3/2 



\\dxP\, 



c f \i {p < 1}P ^- 3/2 d x p\ 2 dx< f 

Jr Jr 



{p<i}' 



d p p(p, p)d x p 



Vpp(p>p)- 



dx 



< C. 



This provides what we desire by way of the following equality: 

ll 1 {p<l}^~~ i '~ 3/2 ^^lli oo (0,T;Lf oc (R)) = C||1{p<1}^P 1 '~- 1 '~ 1/2 ||l~ c (0,T;L 2 (R)) 

< c. 

(4.18) 
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Thus when applying the condition from (2.3) it follows that 77 satisfies 



1 

77 = 017-7+ -. 

Likewise when p > 1 applying (2.5) and (2.8) provides 
d P p(p,p)d x p 



(4.19) 



{p>i}- 



> c\i {p > l} p^'-^ 2 d xP \ = c\i {p > 1} p^- 3 / 2 \\d x p\, 



Vpp{p, m) 

such that integrating over I gives by (4.17) that 



C II 



{p>i}P 



4p>i>- 



Vpp(p,p) u 



dx 



<C K . 



Here this yields 

l|l{ P >i}/^~"^d x p|| L =o (0 , T;i i (/)) = C||l {(0 >i } 9 x p'» 
Thus using the condition from (2.4) establishes 



(7 = 7-0-7- -. 



ll^O.T;!, 1 ^)) < Cif. 

(4.20) 



(4.21) 



In order to complete the proof all that remains is to add (4.18) and (4.20) 
together and apply Minkowski's inequality, which gives 



||^x^(p)||l~(o,T;L 1 (7)) < C K - 



□ 



We are now able to show Proposition 4.1 by applying the preceding results. 
Proof of Proposition 4-1- For t fixed set 



#(1,M(M)) 



for y < 1 
for y > 1 



such that y = p. Next (3.2) together with (2.5) shows that &(y, x) is continuous 
in p uniformly with respect to x, and (4.5) assures that 



/ 



&(p{t,x),x)dx < C. 



Then the hypothesis of lemma 4.2 is satisfied as long as there exists a S > 
such that J^(0, x) > S. But for p < 1 we can apply (2.5) to the form of the 
internal energy (3.2) to see that as p — > we have S{p,p) > C. Likewise when 
p = 1 we see that ${l,p) > C\ for C\ a constant. So we have for a positive 
5 < inf{C, Ci} that the hypothesis of lemma 4.2 is satisfied. Then for any x e M 
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with xq = x from lemma 4.2 there exists anii € I = [x — K , x + K] such that 
p(t,xi) > e. Note that K does not depend on t since § WL ^(p{t 1 x) 1 x)dx does 
not depend on time thanks to (4.5). Then the fundamental theorem provides: 

|l { „<i } p-"(aO| < le^l+^ll^i}^-"!^ 

Since K does not depend on time, lemma 4.4 gives that the right hand side is 
bounded uniformly in x and t. 

For the upper bound, again fix t and now set 

such that y — 1/p. Again (3.2) and (2.5) provide that J?(y,x) is continuous 
in p uniformly with respect to x. Additionally we find that both J?(l/p,x) = 
S{p,p) > C and that ^(0,x) > C\ by applying (2.5) to (3.2), which provides 
an admissible 6. Now, upon defining a function zu = p(l + p)/(p+ 1), then there 
exists a constant C > such that 

g(w,ii)<Cg{p,ii), 

which can be shown using (3.2) and checking the formula for \p — p\ < p < | 
and p > | p thanks to (2.8). Then (4.5) is enough to deduce that 

/ ^(p{t,x)- 1 ,x)dx < C. 

Hence for any x e M we can use lemma 4.2 setting xo — x such that there 
exists an x\ E [x — K, x + K] with p(t, x\) < e _1 . Again notice that K does not 
depend on t since (4.5) is uniform in time. Then by the fundamental theorem 
and lemma 4.4 we obtain 




Again since K does not depend on time, lemma 4.4 gives the right side bounded 
uniformly in x and t which completes the proof of proposition 4.1. □ 

We proceed by showing the important corollary to this proposition. 
Corollary. Assume that (2.3)- (2.5) and (2.8)-(2.9) are satisfied, then 

p&L^^^T-H 1 ^)). 

Proof. Lemma 4.4 provides the appropriate framework. Thus we will show the 
bound separately for the cases p < 1 and p > 1. 
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For p < 1 applying (2.5), (2.8) and (2.9) we calculate 
- v -i ( d x p(p, p) - df,p(p, p)d x p 



= P 



d P p(p, p) 



r,-l ( 9xP(p, P) \ _ 



< Cp-^-^ 2 d xP (p,p) - cp- r >-^ +1 . 

Squaring both sides gives 



(4.22) 



Integrating, applying (2.5) and utilizing Holder's inequality yields, 



I t{ P <i}{d x pfdx < C [ 1 
Jr Jr 



{p<i} 

2 



d x p 



dx 



- 6 U 

\Jr 



{p<i} 



d x p 



Vpp-^ 



< f t {p <_iy P ^-^(?f) dx)\cf 1 {P< _ 1} p^~^(?f) 



< Co l{p<i} 



d x ip(j>) 



VP 



da; -C / l{ p <i} 



d x ip(p) 



VP 



dx 

i 

dx 
2 

dx 



^7(i-«)+l f i {p < 1} \ p -i dx p\ 2 dx) 2 +c^ 1 ^+2 f i^Kp-'d^dx 

Jr J Jr 



<C, 

which concludes the proof for p < 1. 

For the case p > 1 we follow an almost identical calculation, except that now 
after applying (2.5), (2.8) and (2.9); (4.22) becomes 

8 x p° = p°- x d x p 

a -i ( d x p(p, p) - d^p(p, p)d x p 



= P 
<P 



d P p(p, p) 

cr-l ( Cd x p{p, p) 



Cp a '*dp_p(p,p)d x p 



(4.23) 



< c P -^-^ 2 d xP (p, p) - c P °-^ +1 (^) . 

Squaring both sides now gives 

{d xP f < P 2 - 2 ° (cp-^i 2 d x p( P ,p) c P ^-^ (^)) 2 • 
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Again integrating and applying (2.5) with Holder's inequality establishes, 



/ 

JR 



1{ P >i}{d x p) dx <C I l{ P >i} 

2 

. , . - / I) , ■ 1 1 \ 

1 {P>1] 



X Q 



I 

JR 

<C [ 

JR 

!-«) + 5 f t {p > 1} \p-^d xl x\ 2 dx 
JR 

<c, 



d x p 



Vpp° 



dx-C 



( Jr 



P 

l{ P >i} 



dx) + C / 1 



VP 



4p>i} 



dx — Co 



l{ P >i} 

7(1 



d x p 



4~PP ai 



4p>i> 



dx 

! 

2 



Jm. 



which due to Minkowski's inequality completes the proof. □ 
4.2 Bounds for the Velocity 

It is now possible to find bounds on the velocity by applying the uniform bounds 
achieved above. 

Proposition 4.2. Assume that (2.2)-(2.5) and (2.8)-(2.9) are satisfied, then 
u e L 2 (0, T; iJ 2 (R)) and d t u e L 2 (0, T; L 2 (R)). (4.24) 



Proof. First notice that the second estimate in (4.5) in tandem with the uniform 
bounds on the density gives 



IM|l°°(0,T;L 2 (R)) < C. 



(4.25) 



Also notice that the uniform bounds on p applied to (2.2) show that there exists 
a constant C such that v(p, p)^ 1 < C. That is, applying (2.5) and (2.8) to (2.2) 
for p < 1 gives v{p, p) > Cp^^-* 1 so that using the uniform bounds on p provides 

v(p, PY 1 < Cp^-* 1 < Cg^ 1 ^ < C. (4.26) 
For p > 1 it follows in the same way that v(p,p) > Cp* 1 ^^ provides 

v{p, PY 1 < Cp^-" 1 < < C. (4.27) 

Thus for all p we have v{p,p)~ 1 < C, which when applied to (4.5) yields 

l|w||^(0 ; T;/fi(M) < C. (4.28) 

Further, observing the continuity equation with respect to (4.28) implies that 
d t p is bounded in L 2 ((0,T) x R) as denoted in the theorem. We proceed by 
controlling the following form of the momentum equation (after multiplication 
through by p _1 ): 

d t u- d x {p^ 1 v(p,p)d x u) = -ud x u - p^d x p(p, p) - u(p, p)d x ud x p~ 1 . (4.29) 



22 



Compressible Multifmid 



We want to control the right side of (4.29) in such a way as to apply classical 
regularity results for parabolic equations. 

Consider first the second term on the right in (4.29). This term is bounded in 
L°°(0,T; L 2 (R)) as an immediate consequence of proposition 4.1, the corollary, 
and condition (2.5). This follows since (2.5) gives p(p, p) < Cp^ for p > 1 and 
p(p,n) < Cp~' for p < 1. Then we can expand the pressure term as p^~ 2 d x p 
and p^~ 2 d x p, such that for p > 1 the corollary and proposition 4.1 provide that 



/ 

JR 



l{ P >i}|p 1 d x p(p,p)\ 2 dx < C ( esssup \p 21 4 | ( / l {p > 1} \d x p\ 2 dx ) 

\{x£R:p>l} I WR / 



<c, 

and likewise for p < 1 the corollary and proposition 4.1 give 



/ 

Jr 



1{ p <i}Ip 1 d x p(p,p,)\ 2 dx < C ( esssup |p 27 4 | ( / l {p <i } |9 2; /?| 2 dx ) 

\{x(ER:()<l} / \JR J 



< c. 



Minkowski's inequality then provides the result. 

For the third term on the right we again use the fact from above that 
u(p, p)^ 1 < C, and so because of the uniform bounds on p we acquire 

\v(p, p)d x ud x p~ 1 \ < C\d x ud x p\. 

Hence, due to results on parabolic equations (see [58]) we have reduced the prob- 
lem to finding for the third term on the right in (4.29) that p x u x is bounded 
in L 2 (0, T; L 4 / 3 (R)) and similarly for the first term on the right that uu x is 
in £ 2 (0,T;L 4 / 3 (R)). To get this, we adapt a subtle calculation from [71] 
that relies on correctly weighting the norms in order to establish that u x E 
L 2 (0, T; L°°(R)). That is, using Holder's inequality we can write: 

\\ UU x\\l 2 (0 ,T;L 4 / 3 (R)) + IIAe«:e||l 2 (0,T;L 4 / 3 (R)) ^ ^ 

< {II u IIl~(0,T;L 2 (R)) + IIP:c||l~(0,T;L 2 (R))}II u :e|U 2 (0,T;L 4 (R))- 

Now for some function / with constant a £ R we have (f a ) x = af al fx such 
that we may infer by Holder's inequality that 

R(/ 3 / 2 )|| L1(R) < C\\f^ 2 \\ LHR) \\f x \\ Li/3m . (4.31) 
Next we infer a bound in L 8 / 3 (R) given by 

||/ 3/2 IUs/3 (R) < c-ii/^ii^s^ii^a 3 / 2 )!!^, 

which follows since 

WUfV^ > c|l/ 3/2 ll^ (R) . 
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Thus invoking (4.31) we can write 

<C||/II L 2 (R) L 4 / 3 (R) 

< C|l/IU 2 (IR)||/x|li4/3( R ) ) 

where both sides raised to the power n — 2/3 clearly implies that 

II/IIl 4 (r) < cil/lli^li/li^^/a- 

Hence, if we set u x = f then (4.30) leads to 

{ll u IU°°(0,T;L 2 (R)) + IIP:e||l~(0,T;L 2 (R))}I|U:e|U 2 (0,T;L 4 (R)) 

< ^||w||Loo( 0;T ;I, 2 (R))||Wx|li( 3 (o ; T;L 2 (R))ll U ^llL 23 (0,T;VF 1 . 4 /3(R)) 

+ C IIPx||L~(0,T;L 2 (R))l|Wx||£( 3 ( ^ T;i 2( R)) ||Mx||^ 3 (0 ^ T;H/ i,4/3( R) ) 

< C|I u x|Il 2 (0,T;W 1 . 4 / 3 (R))' 

since u and p x are given by (4.28) and the corollary. But then regularity results 
(see theorem 4.2 in Chapter III of [58]) for equations of the form (4.29), given 
the bounds established above and that v{p, p) is a coefficient function satisfying 
uniform parabolicity, imply that since 

1 /3 

\\dxU\\ L 2( T . W l,l/3(M)) < C + C|l U x|| L 2( 1 T;W 1 . 4 /3(R))! 

we have 

ll^a= u llL 2 (0,T;H' 1 . 4 / 3 (R)) < C- (4.32) 

Now, we want to show that 

u x g£ 2 (0,T;L oc (K)). (4.33) 

Indeed for any x € R and t e [0,T] if we set ? = u x from lemma 4.5 (which is 
given following this proof) and notice that 



\\u x {t,x)\\ 2 <2\\u x {t r )\\\ 2{m + \\u xx {t r 



2 



Il 4 / 3 (r) 



for any t G [0,T], then integrating in time gives (4.33). 

It follows as a consequence that the entire right hand side of (4.29) is bounded 
in L 2 (0, T; L 2 (R)). Applying the classical regularity results for parabolic equa- 
tions then yields: 



!|w||l 2 (o,t ; // 2 (r)) < C and ||dtw|U 2 (o,T ; i, 2 (R)) < C. 



□ 



Lemma 4.5. Let <r e i 2 (M) wit/i e ^(R). Then for any x e 

/ r \ 2 

k(x)| 2 <2|M|i 2(R) +2 

w/iere I = [x, x + 1] . 
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Proof. It follows by the fundamental theorem that 

k(x)| < + £ \d x ^\dz < + J \d x c;\dz, 
for any y £ I. Squaring both sides and integrating over R in y yields: 

2 



k(x)| 2 <2|K||| 2(R) +2^\d^\di 



□ 



4.3 Bounds on the Mass Fraction 



All that remains in order to conclude the proof of the existence half of the the- 
orem is to establish the bounds on p,. However, this is now an easy consequence 
of the bounds we have already established above. 

Lemma 4.6. Given proposition 4-1 and 4-2 there exist constants such that, 
IIMk||l°°(0,T;L°°(E)) < C an d l|5tMllL°°(0,T;L 2 (E)) < C. 

Proof. We have from lemma 4.4 that 

||p~ 1/2 0*A*||L~(O,TiL~(R)) < C, (4.34) 

and so thanks to the uniform bounds on the density from Proposition 4.1, this 
yields that d x fi is in L°°(0, T; L°°(R)). Now using (4.25) and the above with 
(1.3) we find that dtfi is in L°°(0, T; L 2 (R)). □ 

4.4 Proof of the Existence Half of the Theorem 

We now apply the preceeding results in §4 in order to prove the existence theo- 
rem. 

Proof of existence half of the theorem. In view of the a priori estimates that we 
have now, the only difficulty that remains is to deal with the fact that v is not 
uniformly bounded by below with respect to p. This is needed to apply the 
short-existence result of Solonnikov (proposition 2.1). To solve this problem let 
us fix any T > 0. Then we define an approximation to v by, 



v(y, z) if y > ^p- 

U{y ' z) = ^{^, z ) if y <s<r> 



where g(T) is defined by proposition 4.1. Now let (p,u,fi) be a strong solution 
of (1.1)-(1.3), where v is replaced by i>\ giving 

d t p + d x {pu) = 0, 
d t (pu) + d x {pu 2 ) + d x p(p,fi) - d x (v{p,jl)d x u) = 0, 
d t {pp) + d x (pup) = 0. 
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By (2.2), (2.5), and (2.8) the approximate function v is bounded from below, 
thus proposition 2.1 provides that such a solution exists for all t € (0,T S ). 
Consider T <T the biggest time such that 

inf05(V))> V- 
x 2, 

Then on [0, T], it follows that v = v. Now assume that T < T. From proposition 
4.1, on [0,f] 

q(T) 

x — z 

which contradicts the fact that T <T. Hence we have constructed the solution 
of (1. 1)-(1 .3) up to time T, and this for any T > 0, which completes the proof. 

□ 



§5 Establishing the Uniqueness Theorem 

Now we address the uniqueness half of the theorem. Thanks to [89] this result 
follows fairly directly. 

Theorem. Let tp"(p), d pp p(p,p) and d p ^p(p 7 p) be locally bounded. Then a 
solution of (l.l)-(l.S) verifying proposition 4-1, proposition 4.2, and lemma 4-6 
is uniquely determined. 

Proof. Let (pi,u\, p\) and (p2, u 2 , P2) be two solutions to the system (1.1)-(1.3), 
and define \ = Mi ~ M2, r = p\ - p 2 , C = "i - Pi = p(pi, Pi) - p(p2, M2) and 
vi = v(pi,pi) — v{p2,p2) such that from (1.1)-(1.3) we can write: 

d t r + d x (piU! - p 2 u 2 ) = 0, 
pid t U! - p 2 dtu 2 + piu x d x ui - p2U 2 d x u 2 + d x pe - d x {vid x u x - v 2 d x u 2 ) = 0, 
dtx + (uid x p! - u 2 d x p 2 ) = 0. 

By rearranging we get 

d t T + d x (TU 1 +p 2 C)^0, (5.1) 

pi(d t C + uid x C, + (d x u 2 ) + T(d t u 2 + u 2 d x u 2 ) 

+ d x p e - d x {v e d x u{) - d x (v 2 d x Q = 0, (5.2) 
d t X + Cd x pi + u 2 d x x = 0. (5.3) 

First let us consider equation (5.1). Here we multiply through by r and 
integrate in x. To begin with, note that the first term on the left satisfies 

f rd t Tdx = \ [ d t T 2 dx. (5.4) 
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For the {tu\) x term we use proposition 4.2 as applied in (4.33) by setting u = u\ 
to see that 

f 1 1 

/ T{T Ul ) x dx < -||T 2 a,;Ui|| L i (K) < 2ll T llz,2(R)ll 5 xWi|| L oo (IR) < Bi(i)l|-r|||2 (R) . 

(5.5) 

For the (p 2 () x term notice that we can write: 

/ Td x (p 2 ()dx < / rp 2 d x (dx + / r(d x p 2 dx 
Jr Jr Jr 

Applying proposition 4.1 and Cauchy's inequality to the first term on the right 
provides, 



rp 2 d x (dx 



<C JjTd x {\dx < C\\T\\ L 2 {m \\( x \\ L2{ 

<C 2 (4 ei )- 1 ||r||i 2(R)+ei ||C ;c ||i 2(R) . 



(5.6) 



For the second term on the right Holder's inequality with the corollary implies 
that 



/ r(d x p 2 dx < ( / |t| 2 (2x) ( / \d x p 2 \ 2 dx) (esssup|C| 
Jr \Jr J \Jr J \ r 



< C 



1/2 



ess sup |C| 



(5.7) 



Now we utilize lemma 4.5 by setting ? — £. Since \(\ < |«i| + \u 2 \ the bounds 
in (4.25) provide that £ e L°°(0, T; L 2 (R)). Furthermore, proposition 4.2 gives 
that since \( x \ 2 < 2\d x m\ 2 + 2\d x u 2 \ 2 we have ( x G L 2 (0, T; L 2 (R)). Thus 
noticing that HCIU 1 ^) < IICIIl 2 (/) since |/| = 1 from lemma 4.5, it follows that 

\{{X)\ < \\C\\L^(R) + ||Cx|Ui(7) < ||CIIl 2 (R) + ||Cr|U»(R), 

allowing us to deduce, 

IICIU~(k) < \\Ch'(R) + II^CIIl^r). 
By Cauchy's inequality this finally yields 

c1MIl2(k)||CIU°°(k) 



< e 2||Cx||i,2 (R) + { — + — 



c 2 c 



Thus combining (5.4), (5.5) and (5.8) allows us to write for (5.1): 



(5.6 



~ I r 2 dx-{e 1 + e 2 } f (d x C) 2 dx 
1 at Jr Jr 



< B 1 (t) 



C 2 

4ei 4e 2 



C 



(» 



L 2 ( 



+ 



IICIIW)) 



(5.9) 
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Next we want to multiply (5.2) through by £ and integrate in R. For the 
first two terms in the first part of (5.2) we find: 

/ PiC($C + «i0*O<&= / ^(dtt 2 +u 1 d x t 2 )dx 
Jr Jr l 

■ [ Pl ( 2 dx- [ ^-(Qtpi + d x (p lUl ))dx (5.10) 

JR JR Z 



ld_ 

Ydt 



For the piC,d x U2 term in (5.2) we use the same calculation given in (5.5) which 
is formulated in (4.33) by setting u = U2 such that, 



/ 

Jr 



piC d x u 2 dx 



< C||Clli 2 ( R )ll^ 2 || L ~ (R) < B 2 (t)||Clli=i(R)- (5-11) 



Now, for the r(dtU2 + Uid x U2) part of (5.2) we utilize a calculation similar to 
that employed for the term in (5.7). Here we simply substitute the d x p 2 term 
from (5.7) with oj = dtU2 + U2d x U2, noting that proposition 4.2 along with (4.33) 
assure that u> is bounded in L 2 (0, T; L 2 (R)). Thus we obtain 

£(i)ll T IU 2 (iR)IICIU~(]R) < e 3||C:i:|||2( R) + e 3 ~ 1 B 3 (t) (||t||| 2 ( R ) + ||CIIl2(r)) > 

(5.12) 

where here B 3 (t) = e 3 B(t)/2 + B(tf/A. 

Next consider the pressure term pi in (5.2). Here set 

/ Qd x pidx = - j \p(pi,m) -p(p 2 ,iJ,2)}d x (dx. 
Jr Jr ' 

The uniform bounds on p along with (2.8) and (2.9) give that \d p p(p,p)\ < C 
and \dnp(p, p)\ < C, and so 

\p(P2,P2)-p(pi,Pi)\<C(\t\ + \ X \). 

Thus 

/ CdxPtdx < C [ (|r| + \ X \)d x (dx, 
Jr Jr 

which gives by Cauchy's inequality, 

/ CdxPtdx < 2e 4 f(d x 2 dx + [ \r\ 2 dx + [ \ X \ 2 dx. (5.13) 
Jr Jr 4e 4 Jr 4e 4 J R 

Finally we consider the viscosity terms in (5.2). For the [yid x u\) x term 

- / C,d x {v t d x u\)dx = j ved x (d x uidx. 
Jr Jr 
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Since ijj"(p), d pp p(p, fi) and d pfl p(p, y) are locally bounded, then from (2.2) we 
have vi < C(\t\ + \x\), which gives 

- / (d x (v e d x ui)dx <C (|r| + \x\)d x C,d x u\dx, 
Jr Jr 



and leads to, 



- / ((vid x ui) x dx < 2e 5 f ( 2 dx + ^— f \d x ui\ 2 T 2 dx + ^— [ \d x u\\ 2 x 2 dx. 
Jr Jr 4e 5 Jr 4e 5 J R 

Next we again use the fact that d x m is bounded in L 2 (0,T; L°°(M)) by (4.33). 
It subsequently follows that, 



-\\C{ved xUl ) x \\ L i {K ) < 2e 5 ||a x C||i,2 (IR) + -^—{\\t\\ L 2 (wl) + ||x|| L2(R) ). (5-14) 

For the (f2Cx)x term we simply multiply through by £ and integrate, yielding 

- / (d x (v 2 d x ()dx = [ v 2 (d x Q 2 dx > C [ (d x Q 2 dx (5.15) 
Jr jr Jr 

when using that v 2 > C. 

Hence combining (5.10)-(5.15) we have: 

~ [ Pl ( 2 dx+(C-e 3 -2e,-2e 5 ) [ \d x (\ 2 dx 

(5 16) 

5 {*<« + ^ + + B ( ||T|li *« + W M ■ 

All that is left is to find a compatible form of equation (5.3). Here we 
multiply through by x and integrate in R such that the first term gives 

/ Xdtxdx = ~ [ X 2 dx. (5.17) 

JR z at JR 

The second term in (5.3) is treated in a similar way as (5.7) and (5.12), where 
here we have 



/ xCdxPidx < f ( \x\ 2 dx\ ( ( \(\ 2 dx\ fesssup \d x m\\ 
Jr \Jr J \Jr / \ r / i- 10 x 



Thus we obtain, 

C 



C||xIU'(r)IICIIl'(R) < ^ (llxll! 2( R) + \\C\\h m ) • (5-19) 
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For the last term in (5.3) we use (4.33) with u = u 2 to see 

Xu 2 d x xdx < C\\x\\ 2 L 2 m \\d x u 2 \\L^(M) < B 5 (t)\\x\\ 2 L 2 {R y 



Thus putting (5.17), (5.19) and (5.20) together yields, 
1 d 



ICII 



Finally, combining (5.9), (5.15) and (5.21) along with defining, 

^ = C - ei + e 2 +e 3 + 2e 4 + 2e 5 , 
#i(t) - Bi(t) + C 2 (4e 1 )" 1 + C 2 (4e 2 )" 1 + C/2, 
# 2 (t) = B 2 {t) + Babies)- 1 + C 2 (4e 4 )" 1 + B 4 (2e 5 )" 1 , 
# 3 (t) = B 5 (t) + C/2, 
*f(t) = 38 1 {t)+38 2 {t) +@s{t) 



yields: 



I A /" + * / i^ci 2 ^ < ^(*) (llxlli a( R) + IICII 

z « r JR JR v 



(5.20) 



(5.21) 



+ h 



LH 



Since proposition 4.1, proposition 4.2 and lemma 4.6 confirm by above that 
s£(t) € L 2 (0,T), and as ^ is positive, then at t = since 

/ ^{t )dx = / r 2 + Pl | t=0 C 2 + Xo^ = 0, 

JR JR 



then Gronwall's lemma gives that J R &{t)dx = over [0,T], which establishes 
that r, £, and \ are eacn zero. □ 
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